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Abstract: A total of 16 Taenia multiceps isolates collected from naturally infected sheep or goats in Gansu Province, Chi- 
na were characterized by sequences of mitochondrial cytochrome c oxidase subunit 1 (coxl) gene. The complete coxl 
gene was amplified for individual T. multiceps isolates by PGR, ligated to pMDIST vector, and sequenced. Sequence 
analysis indicated that out of 1 6 1 multiceps isolates 1 0 unique cox7 gene sequences of 1 ,623 bp were obtained with se- 
quence variation of 0.12-0.68%. The results showed that the coxl gene sequences were highly conserved among the 
examined T. multiceps isolates. However, they were quite different from those of the other 7aen;a species. Phylogenetic 
analysis based on complete coxl gene sequences revealed that T. multiceps isolates were composed of 3 genotypes 
and distinguished from the other 7aen/a species. 
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INTRODUCTION 

Coenurosis, caused by the larval form of Taenia multiceps, is 
a parasitic disease of various livestock species, especially rumi- 
nants. Coenurosis is common in sheep and goats worldwide 
and also has public health concems because of its infection in 
humans [1,2]. T multiceps larval cysts are usually found in ner- 
vous system, including brain and spinal cord of many herbi- 
vores. The presence of cysts typically leads to neurological 
symptoms that cause the infected animals to die after some 
weeks from starvation. However, animals in most cases remain 
normal without clinical symptoms and the disease is diag- 
nosed only after the death of the animals. T multiceps is a tae- 
niid cestode whose adult stage lives in die small intestine of 
definitive hosts such as dogs and other canids [3]. I multiceps 
larvae infection was reported in countries where sheep and 
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goats herding remain important components of local econo- 
mies, including some areas of Europe, Africa, Americas, and 
Asia. In China, T multiceps infection was reported in dogs in 
several provinces with infection rate between 2.3% and 14.6% 
[4,5]. Although coenurosis in animals was reported in more 
than 20 provinces of China, molecular identification of T mul- 
ticeps isolates was rarely recorded. Moreover, the molecular 
identification of T multiceps isolates is of great importance for 
etiological, epidemiological, and phylogenetic studies. 

Mitochondrial (mt) DNA is known to have a faster evolu- 
tionary rate than nuclear DNA and is generally inherited ma- 
temally in almost all metazoans [6,7]. The mt genome is con- 
sidered to be clonal and rarely or never undergoes recombina- 
tion. Sequences generated from the mt genome provide excel- 
lent molecular markers for defining population groups, for 
tracing the genetic history of an individual or a particular 
group of related individuals, and for reconstructing deep- 
branch taxonomic phylogenies. With respect to the molecular 
identification of cestodes, early studies using partial cyto- 
chrome c oxidase subunit 1 gene [coxl) and NADH dehydro- 
genase- 1 gene helped to establish the population and genetic 
relationship of I multiceps and other members of the Taeni- 
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idae [8-10]. However, to our knowledge, there have been rare 
reports on molecular characterization of T. multiceps isolates in 
China. For this reason, we described in this paper an investiga- 
tion of genetic variability among 16 T multiceps isolates from 
Gansu Province, China by using sequences of axl gene as ge- 
netic marker. 

MATERIALS AND METHODS 

Collection of I multiceps isolates 

Between May 2008 and June 2010, a total of 16 isolates of I 
multiceps larval stage were obtained from naturally infected 
sheep or goats in Gansu Province which is situated in the 
northwestem part of China between the northem latitudes of 
32° to 40° and eastem longitudes of 94° to 108°. Sample code, 
geographical origin, and host are shown in detail in Table 1. 
The protoscolices were removed from individual larval cysts, 
washed 3 times with saline and stored at -70°C before extrac- 
tion of genomic DNA. 

Genomic DNA extraction and amplification by PGR 

Genomic DNA was extracted from protoscolices using Ge- 
nomic Tissue DNA Kit (Omega Biotek, Noraoss, Georgia, USA) 
according to die manufacturer's recommendations. Genomic 
DNA concentrations were determined spectrophotometrically 
by GeneQuant 100 and the DNA samples were stored at -20°C 
until used. The primers to amplify the complete coxl gene were 
designed according to the mitochondrial genome sequence of 



Table 1 . Taenia multiceps isolates used in the present study 



Sample code 


Geographical origin 


Host 


Accession number 
of cox7 gene 


JT080526 


Jingtai, Gansu 


Sheep 


JX535567 


JT081008 


Jingtai, Gansu 


Sheep 


JX535568 


JT081204 


Jingtai, Gansu 


Sheep 


JX535569 


JT090318 


Jingtai, Gansu 


Sheep 


JX535569 


JT090104 


Jingtai, Gansu 


Sheep 


JX535570 


PL090603 


Pingliang, Gansu 


Sheep 


JX535570 


JT090115-1 


Jingtai, Gansu 


Sheep 


JX535571 


JT090115-2 


Jingtai, Gansu 


Sheep 


JX535572 


JT090115-3 


Jingtai, Gansu 


Sheep 


JX535573 


JT090603 


Jingtai, Gansu 


Sheep 


JX535573 


JT090331 


Jingtai, Gansu 


Sheep 


JX535574 


JT1 001 24-1 


Jingtai, Gansu 


Sheep 


JX535575 


JT1 001 24-2 


Jingtai, Gansu 


Sheep 


JX535575 


JT1 00202 


Jingtai, Gansu 


Sheep 


JX535575 


YJ1 00610 


Yongjing, Gansu 


Sheep 


JX535576 


ZJCH1 00610 


Zhangjiachuan, Gansu 


Goat 


JX535576 



T. multiceps [11] as follows: PF S'-GACnTCTTTGAGTCATCr- 
TATGTAA-3' and PR 5'-GATrrACAAAATCrACATnTAAA-3'. 
PGR (50 pi) was performed in 10 mM Tris-HCl, 50 mM KCl, 
1.5 mM MgCh, 0.2 mM of each dNTP, 20 pmol of each primer 
and 2 U Taq DNA polymerase (Takara, Tokyo, Japan) under the 
following cycling condirions: after an initial denaturation at 
94°C for 5 min, then 94°C for 1 min (denaturation), 55°C for 1 
min (annealing), 72°C for Imin (extension) for 35 cycles, fol- 
lowed by a final extension at 72°C for 10 min. Samples without 
genomic DNA were included in each amplification mn as nega- 
tive controls. PGR produrts were separated on a 1% agarose gel 
and detected by ethidium bromide staining. 

DNA sequencing and analysis 

All of the amplicons were ligated into pMD-18 T vector (Ta- 
kara) and transformed to Escherichia coli DH5a competent 
cells. The positive colonies were sequenced with ABI PIUSM 
3730 Genetic Analyzer, and sequences were analyzed with 
DNAtools software. Further comparison with taeniid coxl gene 
sequences available in GenBank was made using the BLAST 
and Clustal W software. Phylogenetic analyses were performed 
using the neighbor joining (NJ), maximum parsimony (MP), 
and minimum evolution (ME) methods, and the Mega 5.05 
program. The bootstrap value was set as 1,000 replications 
with a cutoff value of 50%. 

RESULTS 

PGR amplification of cox7 gene 

In the present study, 1 to 3 larval cysts were observed in the 
brain of naturally infected sheep, while 7 cysts were found in 
the brain of the infected goat. The characteristics and mor- 
phology of T. multiceps larvae were seen in all of the cysts con- 
taining dear fluid with numerous fertile protoscolices (Fig. 1). 
Using 16 isolates of I multiceps protoscolices DNA as tem- 
plates, a band of approximately 1.7 kb fragment was generated 
after PGR amplification of coxl genes. Each of the PGR prod- 
ucts was ligated to pMD-18 T vector, and 3 colonies for each 
sample were sequenced on both direction and tiie consensus 
sequences were obtained. 

Sequence analysis of cox1 gene 

After removal of the primer sequences, the coxl gene from 
each isolate was 1,623 bp in length encoding 540 amino acids. 
Sequencing results indicated that out of 16 T multiceps isolates 
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10 unique coxl gene sequences were obtained with sequence 
variation of 0.12-0.68% and a total of 19 nucleotide variation 
sites. Some of the isolates had identical coxl gene sequences, 
such as isolates IT081204 and )T090318, isolates ]T090104 and 
PL090603, isolates JT09115-3 and JT090603, isolates YI100610 
and ZJCH100610, and isolates JT100124-1, 11100124-2 and 
JT100202. The G+C contents of the coxl gene sequences of I 
multiceps isolates were 29.9-30.4%. The percent identiries 
among T. multiceps isolates from Gansu Province were more 
than 99% and showed a lesser degree of var iability. Jingtai iso- 
lates JT081204 and 1T090318 were idemical with the corre- 
sponding sequence available in GenBank (GQ228818). When 
the 10 unique coxl gene sequences were translated into amino 
acids sequence, 7 unique protein sequences were obtained 
with sequence variation of 0.19-0.56% and a total of 6 amino 
acid variation sites. Amino acid sequences of isolates JT080526, 
JT090104, JT090115-2, and YI100610 were identical to the cor- 
responding sequence available in GenBank (ACK55903) and 
that of isolate JT081204 was identical to I multiceps coxl avail- 
able in GenBank (ACS37254). 

When Blast analysis with coxl gene sequences was performed, 
more than 20 T multiceps partial coxl gene sequences from 346 
to 445 bp were also available in GenBank. Partial coxl gene se- 
quences from the 10 unique sequences in this study could be 
further grouped into 4 unique sequences with variation rate of 
0.25-0.75%. The first group included JT080526, )T090104, and 
JT100124-1, and was identical to the Turkish T. multiceps isolate 
tmtrOl (EF393620). The second group consisted of ]T081204, 
jT090115-3, and JT090331, and was identical to die Italian T. 




multiceps isolate Tm2 (DQ309768). The third group included 
JT081008, 1T090115-2, and YI100610, and tiie last group was 
JT090115-1. Gomparison of all the partial coxl gene sequences 
revealed 10 groups with variation rate of 0.25-4.55%. 

Phylogenetic analysis of taeniid cox1 gene 

The complete coxl gene sequences from the other Taenia mt 
genomes available in GenBank were selected and aligned. Ac- 
cordingly, a phylogenetic tree was constructed with coxl gene se- 
quences from this study and previously published mt sequences 
of related Taenia species, including I multiceps (GQ228818 and 
FJ495086), T saginata (AY684274), T asiatica (AF445798), T 
solium (AB086256), T hydatigem (GQ228819 and FJ518620), I 
crassiceps (AF216699), T pisiformis (GU569096 and NG013844), 
and I taeniaefomis (JQ663994 and FJ597547), witii E. granulo- 
sus (AF297617) as the outgroup. The coxl gene sequences were 
1,623 bp in lengtii for all I multiceps isolates, but it was 1,620 
bp for T saginata, T. asiatica, T. solium, T. hydatigena, and I pisi- 
formis. However, it was 1,614 bp and 1,635 bp for T crassiceps 
and T. taeniaeformis, respectively. 

Phylogenetic reconstmction of coxl gene from Taenia spedes 
using the NJ, MP, and ME analysis showed similar topology; 3 
genotypes were revealed in the examined T multiceps isolates 



■ ♦T.multiceps JX635673 

♦ T.multiceps JX535574 

♦ T.multiceps JX535571 
T.multiceps GQ228818 
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Fig. 1 . T. multiceps larval cyst from the brain of a naturally infect- 
ed sheep. 



Fig. 2. Phylogenetic relationship among the examined Taenia 
species inferred by neighbor joining (NJ), maximum-parsimony 
(IVIP), and minimum-evolution (ME) analyses based on mitochon- 
drial cox7 sequences, using Ecliinococcus granulosus as the 
outgroup. The numbers along branches indicate bootstrap values 
resulting from different analyses in the order: NJ/MP/ME. Values 
less than 50 are given as 
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(Fig. 2). The first genotype consisted of )T090115-3, JT090115- 
1, 1 multiceps GQ228818, JT081204, and ]T090331. The sec- 
ond genotype included JT081008, JT090115-2, and YJ100610. 
The third genotype comprised of JT090104, JT101124-1, 
JT080526, and T multiceps FJ495086. The axl gene sequences 
of T multiceps were more similar to those of T. saginata and I 
asiatica, suggesting that the coxl gene could be a good genetic 
marker for molecular identification of Taenia species. 

DISCUSSION 

Cestodes of the family Taeniidae (taeniids) occur as adult 
tapeworms in the small intestine of carnivorous definitive 
hosts, and are transmitted to specific intermediate mammalian 
hosts where they develop as fluid-filled metacestodes in tis- 
sues. A number of Taenia species including T multiceps are of 
medical and/or economic significance [12]. For the molecular 
characterization and phylogeny study of taeniid cestodes, 
many genetic markers have been used, including 28S rDNA, 
coxl, nadl, md4, ITS rDNA, and nuclear- protein coding genes 
such as rphz, pepck, and pold [8-10,13-18]. The coxl gene is one 
of the most used markers for Taenia identification. Varcasia et 
al. [19] analyzed genetic variation within 40 T multiceps iso- 
lates obtained from various locations of Sardinia. Partial se- 
quences of nadl and coxl showed differences ranging from 
1.27 to 2.54% and from 0.22 to 0.67%, respectively, and a 
conclusion was made that Sardinian sheep samples had at 
least 3 specific genetic variants [19[. However, comparison be- 
tween the partial coxl sequences of the T. multiceps isolates 
from infected cattie from Erzumm Province, Turkey and other 
T multiceps isolates available in GenBank showed nucleotide 
differences ranging from 0.2 to 2.6% [20]. 

Although coenurosis in sheep and goats were reported in 
China, the molecular identification of T. multiceps was rarely 
conducted. Recentiy, the complete mt genome of I multiceps 
from Chinese isolates was reported and provided usefijl mark- 
ers for studying the systematics, population genetics, and mo- 
lecular epidemiology of Taenia ] 13,21]. Partial coxl genese of 3 
larval isolates from goats in Hunan Province were sequenced 
and showed 97.9% identity to the T. multiceps coxl gene 
(EF393620), indicating that these parasites in goats belong to 
the larvae of T. multiceps [22]. The partial coxl gene sequences 
in this study could be divided into 3 types with variation rates 
of 0.25 to 0.75%, which is similar to the previous report 119]. 
The phylogenetic tree based on coxl gene sequences showed 



that all T multiceps isolates constitute a branch consisting of 3 
genotypes. T multiceps, T. saginata, and T. asiatica showed the 
closest relationships, followed by I solium, T. pisiformis, T. hy- 
datigena, T. crassiceps, and the otiier Taenia species. Therefore, 
the coxl gene sequence can be used as a molecular marker for 
the interspecific identification oi Taenia species. 

When partial sequences of the coxl gene between I asiatica 
and T. saginata were compared, the sequence differences varied 
from 2% in the conserved region (nucleotide position of 742- 
1,107) to 7% in the variable region (positions 25-340) ]23]. It 
was suggested that partial coxl sequences might give a bias to 
analysis, but the complete sequence data should provide a 
more reliable resuh. In this study I multiceps metacestode (lar- 
val form) from naturally infected sheep or goat in 4 geographi- 
cal origins in Gansu Province were genetically characterized 
for the first time with complete coxl gene sequences. We have 
detected some polymorphisms within the common 'sheep' 
isolates. It is evident that intraspecific variation existed in I 
multiceps, whether such minor genetic variation has biological 
or epidemiologic significance warTants fiarther study In addi- 
tion, for studying the population genetics of I multiceps from 
sheep and goats in China, more isolates should be collected 
and analyzed. Further investigations will require more com- 
plete coxl gene sequences from I multiceps isolates from dif- 
ferent endemic regions in China and will determine whether 
there are morphological, developmental, or pathological dif- 
ferences between the genotypes. 
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